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I .  INTRODUCTION 


Carbon-carbon  (C-C)  composites  have  established  themselves  as  impor¬ 
tant  high-temperature  materials  for  spacecraft  applications.  Proposed  new 
aerospace  designs  such  as  those  in  the  National  Aerospace  Plane  (NASP)  and 
in  advanced  turbine  engines  also  utilize  C-C  composites  extensively. 

Some  C-C  composites  are  graphitic  materials.  Such  materials,  includ¬ 
ing  carbon-fiber-reinforced  carbon-matrix  composites,  are  generally  manu¬ 
factured  at  temperatures  exceeding  2400° C.  Glassy-carbon-matrix  C-C 
composites  also  are  frequently  heat-treated  to  relatively  high  tempera¬ 
tures,  above  1700°C,  to  achieve  their  maximum  mechanical  properties.1  When 
such  advanced  composite  materials  are  used  for  applications  requiring  high 
thermal  stability,  thermal  strain  becomes  an  important  engineering 
property . 

The  bulk  thermal  strain  of  C-C  composites  is  ideally  described  in 
terms  of  the  properties  of  the  constituent  materials.  Improved  composites 
can  then  be  designed  without  the  laborious  trial -and-error  of  manufacturing 
test  composites.  However,  C-C  composites  are  complex,  and  this  type  of 
predictive  design  is  difficult  to  apply.  Thermal  strain,  the  reversible 
dimensional  change  associated  with  heating,  is  an  approximately  linear 

p 

expansion  for  most  solids  but  not  for  graphite.  Graphite's  basal  plane 
has  a  minimum  strain  near  400°C,  and  there  is  a  large  crystallographic 
thermal  expansion  anisotropy  between  the  basal  plane  and  the  c-axis.  These 
factors  make  precise  thermoelastic  modeling  of  a  composite  difficult. 

Also,  data  are  lacking  on  the  pyrolysis  and  high-temperature  behavior  of 
precursors  and  matrix  carbons.  During  the  manufacture  of  a  C-C  composite, 
the  carbon  matrix  precursor  is  pyrolyzed  and  carbonized  to  convert  it  to 
carbon  or  graphite.  In  the  process,  the  matrix  densifies  and  a  complex 
crack  network  and  stress  field  form  in  the  matrix  and  at  the  matrix- 
filament  interfaces, ^  making  it  difficult  to  form  an  a  priori  estimate  of 
the  final  thermal  strain  of  a  given  composite  design.  Further  investiga¬ 
tion  of  the  matrix  pyrolysis  behavior  and  its  correlation  with  bulk  thermal 
strain  is  therefore  required. 
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A  problem  arises  in  laboratory  investigations  of  such  phenomena  over 
the  2000-to-3000°C  range:  The  upper  temperature  limit  of  most  coEBercially 
available  high -temperature  dilatoroeters  is  only  about  2000°C.  We  therefore 
constructed  and  tested  a  high-temperature  graphite  dilatometer  that  can 
accurately  measure  thermal  strain  phenomena  at  temperatures  as  high  as 
3000°C.  Before  this  apparatus  was  constructed,  high-temperature  dilatomet- 
ric  and  rheological  measurements  could  not  be  carried  out  at  The  Aerospace 
Corporation,  and  such  facil'ties  are  generally  not  available  at  commercial 
testing  laboratories. 

This  report  describes  the  dilatometer  and  is  organized  as  follows. 

The  primary  components  of  the  dilatometer  system  are  described  in  Section  II. 
The  theoretical  and  practical  aspects  of  dilatometer  calibration  and  the 
method  of  data  reduction  are  discussed  in  Section  III.  Finally,  initial 
experiments  to  measure  thermal  strain  on  several  carbonaceous  materials 
with  the  dilatometer  are  summarized  in  Section  IV. 
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II.  DESCRIPTION  OF  APPARATUS  AND  SAFETY  CONSIDERATIONS 


The  dilatometer ,  diagrammed  in  Fig.  1,  consists  of  a  sample  tube* push- 
rod  and  a  cage  assembly  that  are  fabricated  from  a  single  billet  of  poly- 
crystalline  graphite  to  ensure  uniformity  of  thermal  expansion.  The  sample 
cage  can  be  removed  from  the  sample  tube  so  that  other  cages  can  be  fabri¬ 
cated  and  installed  to  accommodate  samples  of  different  sizes.  The  sample 
pushrod  extends  coaxially  through  the  sample  tube  to  make  contact  with  the 
sample.  The  mass  of  the  pushrod  (-10  g)  serves  to  keep  it  in  contact  with 
the  sample  during  the  test.  A  linear  variable  differential  transformer 
(LVDT),  whose  core  is  fixed  to  the  pushrod,  measures  the  displacement  of 
the  pushrod  (as  the  sample  expands  upon  heating)  and  sends  an  electrical 
signal  to  a  data  acquisition  system.  The  heating  elements,  which  are  also 
made  from  graphite,  are  located  coaxially  around  the  sample  cage. 

The  furnace  (manufactured  by  Astro  Division  of  Thermal  Technologies 
Inc.)  consists,  from  the  heating  elements  outward,  of  a  graphite  retort 
around  which  is  packed  thermal  insulation  (lampblack  and  graphite  felt), 
encased  in  an  airtight,  water-cooled  aluminum  housing.  The  interior  of  the 
aluminum  housing  is  purged  with  an  inert  atmosphere  during  the  test  to 
prevent  oxidation  of  the  insulation  and  retort,  heating  elements,  sample 
holder  assembly,  and  sample.  The  furnace  case  and  the  LVDT  are  cooled  by 
flowing  water  from  the  laboratory  cooling-tower  water  supply. 

The  cooling  water  has  an  emergency  city-water  (tap  water)  backup.  If 
the  cooling-tower  water  flow  is  interrupted,  e.g.,  by  a  power  outage,  two 
valves  can  be  manually  opened  to  provide  city-water  cooling  to  the  furnace 
case.  The  cooling-tower  water  lines  are  monitored  by  a  flow  sensor  that 
sounds  an  alarm  (100-dB  intensity)  and  shuts  off  the  furnace  power  when  the 
cooling-water  flow  rate  drops  below  normal.  Following  such  an  interrup¬ 
tion,  the  time  allowed  at  a  given  temperature  before  the  aluminum  furnace 
housing  fails  from  overheating  has  been  estimated  by  heat  flow  calculations 
(see  Fig.  2).  At  a  maximum  operating  temperature  of  2500t>C,  5  min  would  be 


7 
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Fig.  2.  (a)  Heat  flow  out  of  furnace  case  as  measured  from  change  in 

the  temperature  of  the  cooling  water  flowing  at  a  known  rate, 

Q  =  me  AT:  m  is  the  coding-water  flow  rate,  c  is  the  heat 
capacity  of  water,  AT  is  the  temperature  rise  of  the  cooling 
wated  and  Q  is  the  heat  flow  into  the  furnace  shell  from  the 
interior,  (b)  Approximate  times  for  the  exterior  of  the 
furnace  to  reach  melting  temperature.  Note  also  that  these 
calculations  assume  constant  heat-flow  rate,  i.e.,  that  the 
interior  of  the  furnace  stays  at  constant  temperature  after 
water  interruption,  which  it  does  not.  The  actual  time  to 
failure  is,  however,  regarded  as  50%  shorter  than  indicated  in 
order  to  allow  a  safety  margin. 
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available  foi  the  operator  to  shut  down  the  power  so  that  no  potential 
safety  problem  is  presented.  The  possibility  of  a  cooling-water  failure 
dictates  that  furnace  operation  be  carried  out  only  during  the  normal 
working  shift,  limiting  run  length  to  about  8  hr.  A  "tag-team"  system  is 
used  to  ensure  the  continuous  presence  of  a  qualified  indi vidua-  in  the 
laboratory  during  that  shift.  Figure  3  is  a  diagram  of  the  cooling  system 
with  the  emergency  bypass  system  noted. 

Sample  dilation  and  temperature  data  are  recorded  simultaneously  by  a 
multitasking  laboratory  computer,  a  Hewlett-Packard  Integral  PC.  The  data 
from  completed  runs  are  sent  as  ASCII  files  for  analysis  to  a  VAX  11  785 
mini -mainframe  computer  via  the  Intelligent  Data  Switch  (IDS).  The  Inter¬ 
active  Data  Language  (IDL)  programing  language,  used  for  data  analysis  on 
the  VAX  computer,  greatly  reduces  programming  and  data  analysis  time  while 
retaining  the  inherent  flexibility  of  user-designed  data  analysis  software. 
The  BASIC  language  code  for  data  acquisition  is  in.luded  in  Appendix  A,  and 
the  IDL  code  for  data  analysis  is  in  Appendix  B. 

The  temperature  of  the  furnace  and  that  of  the  test  sample  are 
measured  by  two  different  systems:  For  temperatures  below  1000°C,  a 
Chromel-Alumel  thermocouple  is  positioned  in  the  sample  cage,  just  below 
the  sample.  For  the  900-to-3000°C  range,  optical  pyrometry  is  used  and 
requires  two  separate  pyrometers:  a  total -emission-type  pyrometer,  which 
covers  the  range  1000  to  1550®C  and  measures  the  total  flux  of  radiation 
from  a  specific  area  of  the  sample;  and  a  so-called  two-color  pyrometer, 
which  reads  from  1550  to  3000°C  and  measures  the  intensity  ratio  of  two 
adjacent  near- infrared  wavelengths.  For  both  pyrometers,  the  sample's 
emissivity  is  assumed  to  approximate  that  of  a  blaekbody.  For  maximum 
accuracy,  a  small  cavity-radiator  serves  as  a  target  for  the  pyrometer 
measurements.  The  cavity-radiator  was  used  to  reduce  the  undesirable 
reflectivity  that  the  system  may  have  on  the  sample  so  that  better  accuracy 
can  be  achieved.  One  side  of  the  radiator  is  located  within  1  mm  of  the 
sample  during  the  run  to  maintain  near-isothermal  conditions.  As  the 
furnace  temperature  passes  1000°C  during  a  typical  run,  the  thermocouple  is 
withdrawn  manually.  The  temperature -measurement  function  automatically 
passes  from  the  thermocouple  to  the  optical  pyrometers  at  about  950°C. 
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Fig.  3.  Diagram  of  the  cool ing- tower  water  system  for  the 

dilatometer  furnace.  If  the  system  fails,  valves  A  and 
B  are  closed  and  valves  D  and  C  are  opened  (in  that 
order)  to  provide  a  city-water  backup. 
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The  selection  of  the  particular  inert  purge  gas,  helium,  for  the 
diiatometer  is  important  for  very-high-tempe. ature  applications  with  an 
electric  field  present.  At  temperatures  above  about  25^0*0,  argon,  the 
most  commonly  used  of  the  noble  gases,  begins  to  ioni2e  and  conduct 
electricity.  The  location  of  both  electrodes  at  the  same  end  of  the 
heating  element  in  this  particular  furnace  design  results  in  a  large 
electric  potential  gradient  there  (-50  V/cm}.  When  attempts  were  made  to 
heat  the  furnace  above  in  an  Ar  atmosphere,  a  spontaneous  short 

developed  across  the  heating  elements  and  the  furnace  would  cool.  The 
short  carried  about  300  A  at  10  V  and  was  attributed  to  thermally  assisted 
ionization  of  the  Ar.  When  helium  was  substituted  in  the  furnace 
atmosphere,  the  short  disappeared  completely. 

The  useful  upper  temperature  range  of  the  diiatometer  is  dictated 
primarily  by  the  lifetime  of  the  polycrystalline  graphite  heating  element. 
At  3000°C,  the  lifetime  of  the  element  is  less  than  30  min;  below  this 
temperature,  the  life  of  the  element  increases  substantial ly .  Wachi  and 
Gilmartin^  found  that  at  high  temperatures  (>2900°C),  the  carbon  binder 
phase  in  a  polycrystalline  graphit;  begins  to  evaporate  preferentially, 
generating  gaps  between  the  graphite  grains;  the  graphite  grains  are  then 
(apparently)  ejected  from  the  degrading  material.  Electric  arcs  that  may 
form  in  the  resulting  gaps  in  the  diiatometer  heating-element  material 
might  also  accelerate  the  erosion  process.  The  service  life  of  the  heating 
element  when  the  furnace  is  run  to  3000°C  is  about  three  weeks,  and  that  of 
the  furnace  retort  is  about  three  months.  Upper  temperatures  are  therefore 
usually  limited  to  2850°C. 
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III.  CALIBRATION  OF  THE  INSTRUMENT 


A.  THEORETICAL 

The  accuracy  of  the  single-pushrod  dilatometer  design  is  limited  by 
the  thermal  expansion  of  the  materials  out  of  which  the  dilatometer  is 
fabricated.  Vitreous  silica  is  the  preferred  material  for  high-precision 
dilatometers  because,  for  many  applications,  its  coefficient  of  thermal 
expansion  (CTE)  contributes  negligibly  to  the  overall  expansion  of  the 
sampie/dilatometer  system.  However,  vitreous  silica  dilatometers  must  be 
corrected  for  high-accuracy  measurements  on  materials  whose  thermal  expan- 
sions  are  close  to  that  of  silica  glass  itself  {about  5  *  10"'  °C  ). 

Silica  cannot  be  used  in  this  dilatometer,  however,  because  the  useful 
maximum  temperature  for  silica  is  only  about  1100°C.  The  dilatometer 
discussed  here,  therefore,  is  constructed  of  a  fine-grained  polycrystall ine 
graphite  (SPEER,  Inc.)  to  permit  high  temperatures  to  be  investigated.  The 
average  CTE  of  such  graphites  ranges  from  about  3  to  7  «  10"^  “C"1  from 
room  temperature  (RT)  to  3000°C,  which  is  about  an  order  of  magnitude 
greater  than  that  of  silica.  Consequently,  the  dilatometer  must  be  care¬ 
fully  calibrated  to  remove  the  apparent  CTE  contribution  of  the  probes  in 
the  sample  pushrod/tube  assembly. 

The  procedure  we  used  for  correcting  the  dilatometer  data  arises  from 
the  mechanics  of  the  dilatometer  system  and  is  generalized  for  single- 
pushrod  dilatometers;  it  will  now  be  described.  Thermal  strain  occurs  at 
specific  locations  in  the  dilatometer,  as  indicated  in  Fig.  4a.  For 
simplicity,  the  figure  includes  only  the  primary  structural  parts  of  the 
dilatometer.  The  three-point-bend  fixture,  which  replaces  the  dilation 
sample  in  the  hot  zone  for  flexural  studies,  is  shown  in  Fig.  4b.  The  same 
correction  procedure  can  be  used  in  the  flexure  mode  if  an  inflexible 
standard  sample,  such  as  pyrolytic  graphite,  is  used  for  the  calibration. 

The  sample  and  the  part  of  the  pushrod/tube  assembly  that  are  in  the 
hot  zone  are  heated  during  a  test  and  undergo  thermal  strain.  The  total 
strain  is  measured  by  the  LVDT  at  C.  The  contribution  of  the  "overlapping" 
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(a)  Schematic  of  single-pushrod  dilatooeter. 

(b)  Detail  of  three-point-bend  fixture.  (F  is  force.) 
During  use,  the  fixture  is  positioned  in  sample  zone  A 
of  (a). 


portion  of  the  probes  (in  region  B)  to  the  total  measured  strain  is  2ero  if 
the  materials  that  compose  them  are  identical  in  thermal  expansion,  and  if 
the  temperature  gradients  along  the  tube  and  the  pushrod  are  the  same. 
However,  the  thermal  strain  of  the  tube  and  pushrod  do  not  cancel  in  the 
region  where  the  sample  is  located  (D);  the  thermal  strains  from  both  the 
sample  (A)  and  the  sample  tube  in  region  D  contribute  simultaneously  to  the 
net  thermal  strain  measured  at  the  LVDT.  Therefore,  the  LVDT  actually 
measures  the  thermal  strain  difference  between  the  sample  length,  A,  and 
the  corresponding  length  of  the  sample  tube,  D.  The  arrangement  of  the 
probes  requires  that  the  thermal  strain  of  the  length  of  the  sample  tube  in 
region  D  be  added  algebraically  to  the  total  measured  thermal  strain, 
because  the  strain  of  this  part  of  the  tube  acts  to  move  the  center  of  mass 
of  the  sample  downward,  away  from  the  LVDT.  The  thermal  strain  character¬ 
istic  of  the  sample  tube,  at  least  in  the  region  overlapped  by  the  sample, 
must  first  be  precisely  known  before  this  can  be  done. 

The  characteristic  uniform  thermal  strain  of  the  sample  tube  should  be 
ascertainable  by  measuring  the  thermal  strain  of  a  well-characterized 
reference  material.  However,  it  was  found  that  the  thermal  strain  of  the 
graphite  sample  tube  could  not  be  normalized  with  respect  to  change  in  the 
sample  length  (see  Fig.  7).  Therefore,  measurements  of  a  single  reference 
sample  of  fixed  length  were  not  adequate  to  characterize  the  tube  for  all 
sample  lengths.  (The  correction  for  this  effect  will  be  discussed  in  the 
next  section.)  The  difficulty  appears  to  be  due  to  the  large  and  unequal 
thermal  gradient  in  the  sample  tube  and  pushrod  at  the  edge  of  the  hot 
zone,  and  the  fact  that  different  lengths  of  the  pushrod  enter  the  hot  zone 
for  different  sample  lengths. 

To  bypass  this  problem  without  extensively  characterizing  the  baseline 
associated  with  various  sample  lengths,  the  sample  length  was  standardized 
at  1  cm.  The  measured  strain  at  the  LVDT,  then,  can  be  thought  of  as  the 
sample  strain  plus  the  dilatometer 's  instrumental  response: 

NET  ,n 

“SAMPLE  "  “SAMPLE  *  “MACHINE  '  ’ 
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which  can  be  rewritten 


v  . 

“SAMPLE  =  '“SAMPLE  "  °STD;  *  °STD  *  “MACHINE  ' 

where  aSTp  is  the  known  expansion  of  a  standard  reference  material.  A 
sample  of  the  standard  reference  material  can  be  run  in  the  dilatometer, 
yielding,  by  Eq.  (1), 


°STD  =  “STD  +  “MACHINE 
NET 

The  net  strain  of  the  standard  can  be  subtracted  from  that  of  the 

unknown,  Eq.  (2),  to  yield 


NET 

“CORRECTED  s  “SAMPLE  "  “STD 


(4) 


Since  aSTD  was  assumed  to  be  known,  it  can  be  added  (as  previously 
mentioned)  to  this  result  to  yield  the  true  thermal  strain  of  the  unknown 
sample. 

The  material  for  a  standard  sample  should  be  selected  because  it  is 
either  precisely  characterized  or  it  possesses  a  low  enough  CTE  that  any 
imprecision  in  the  characterization  of  its  CTE  is  of  little  significance. 
Practical  candidates  are  silica  (fused  quartz  glass),  Zerodur-type  glass 
ceramics  (5  *  10"®  °C“1),  and  titanium  silicates  (<  5  *  10“^  °C-1)- 
Partially  crystallized  glasses,  such  as  Zerodur,  are  limited  to  about 
600°C.  Silica  can  be  subjected  to  almost  1100°C  and  can  be  made  more 
refractory  and  lower  in  CTE  by  the  addition  of  up  to  10J  titania  (ULE).® 
The  melting  ranges  of  titanium  silicate  glasses  have  not  been  clearly 
identified  in  our  laboratory,  but  such  glasses  appear  to  withstand 
temperatures  of  about  1400°C  without  melting,  although  devitrification 
occurs  at  the  sample  surface.  In  general,  the  above  materials  offer  such 
low  expansions  that  even  an  error  as  large  as  50>  in  their  CTE  data  repre¬ 
sents  an  order-of-magnitude-lower  error  than  is  expected  from  the  graphite 
dilatometer.  (Silica- based,  single-pushrod  d i la tome ter s  of  this  type  are 
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practically  limited  to  measurements  of  the  order  of  10"^  m  due  to 
mechanical  errors  in  the  apparatus.^)  National  Bureau  of  Standards  (NBS, 
now  referred  to  as  NIST,  National  Institute  of  Standards  and  Technology) 
standard  reference  material  (SRM)  738  (sapphire),  which  is  no  longer 
available  from  NIST,  can  also  be  used  since  it  is  very  accurately  cali¬ 
brated,  but  its  temperature  range  is  limited  to  about  1550°C. 

B.  PRACTICAL  ASPECTS  OF  D1LAT0METER  CALIBRATION 

Polycrystalline  graphite  has  a  very  high  use  temperature,  and  conse¬ 
quently  it  has  been  employed  as  a  standard  material,  in  other  work,  after 
careful  optical  calibration.^  However,  pyrolytic  graphite  (PG)  has  a  lower 
CTE  that  is  considered  to  make  it  a  better  standard  material  because  errors 
in  estimating  its  CTE  would  be  correspondingly  low.  The  theoretical 
(crystallographic)  expansion  of  graphite  parallel  to  the  basal  plane  is 
about  that  of  silica  (1  *  10"^  °C_1),  a  fact  that  favors  its  selection; 
however,  owing  to  the  crystalline  imperfections  inherent  in  PG,  an  expan¬ 
sion  about  two  times  higher  than  this  is  usually  observed. 

For  the  above  reasons,  a  standard  PG  sample  was  used  in  this  dilatom- 
eter  and  was  developed  in-house  as  follows.  A  7-cm-long  sample  of  PG  was 
pre-annealed  at  2750°C  for  2  hr  in  a  graphite  tube  furnace,  then  reheated; 
its  thermal  strain  was  measured  optically  (above  850°C)  using  traveling 
telemicroscopes  that  were  sighted  on  fiducial  marks  in  the  ends  of  the 
rod.  A  1-cm  piece  was  excised  from  this  sample  for  use  as  a  standard.  It 
was  then  calibrated  to  750°C  in  a  DuPont  9^3  thermomechanical  analyzer 
(silica  dilatometer) ,  which  was  referenced  to  NBS  silica  (SRM  739).  These 
thermal  strain  data  were  then  consolidated  to  yield  the  thermal  strain 
behavior  of  the  sample  over  the  entire  temperature  range  of  the  dilatom¬ 
eter.  A  polynomial  curve  was  then  fit  to  the  data  for  use  in  subsequent 
data  analysis.  Figure  5  plots  the  data  from  this  standard  and  the  polynom¬ 
ial  fit.  Published  data  for  the  thermal  strain  of  ATJ-S  polycrystalline 

O 

graphite  are  plotted  in  curve  (b),  Fig.  5,  for  reference. 
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Fig.  5.  Thermal  strain  data  from  (a)  optical  calibration  of  the  in-house- 
developed  pyrolytic  graphite  standard,  (b)  published  data  for 
ATJ-S  polycrystalline  graphite,0  and  (c)  theoretical  basal-plane 
expansion  of  graphite.  The  dashed  curve  (d)  is  the  corrected  PC 
data  from  the  in-house  standard. 
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Measurements  of  the  NBS  738  sapphire  sample  in  the  dilatometer  after 
this  calibration  procedure  indicated  the  need  for  applying  the  following 
correction  to  the  optical  calibration  data  for  the  PG  standard.  The  CTE 
data  for  a  sample  of  NBS  738  sapphire,  together  with  the  published  data,^ 
are  plotted  in  Fig.  6.  None  of  the  materials  involved,  in  either  the  stan¬ 
dards  or  the  equipment,  undergoes  phase  changes  during  heating,  and  all  are 
considered  homogeneous  and  isotropic.  Therefore,  it  was  assumed  that  these 
materials  display  expansion  characteristics  described  by  a  function  that  is 
continuous,  strictly  increasing,  and  has  first  and  second  derivatives  that 
are  also  continuous  and  strictly  increasing  (i.e.,  the  curve  is  very 
"smooth,"  having  no  inflections).  This  is  a  good  assumption  since  graphite 
is  nonpolymorphic  at  the  temperatures  and  pressure  used,  and  pre-heat 
treatment  to  high  temperatures  renders  it  stable  at  lower  tempera¬ 
tures.1®’11  In  order  to  make  these  smoothness  assumptions,  the  dilatometer 
had  been  preconditioned  by  heating  to  3200°C  for  0.5  hr  prior  to 
calibration. 

The  deviation  of  the  measured  sapphire  standard  data  from  the 
published  data  increased  very  smoothly  (see  deviation  data  in  Fig.  6a), 
thereby  indicating  that  a  linear  error  was  present  in  the  optical  calibra¬ 
tion  data  of  the  PG  standard.  The  indicated  correction  for  this  deviation 
was  made  to  the  optical  calibration  data  by  extrapolating  the  deviation  to 
3000°C  and  subtracting  the  resultant  curve  from  the  PG-standard  data 
(Figs.  5  and  6) . 

To  observe  the  dilatometer 's  fine-scale  instrumental  response,  samples 
of  ATJ-S  polycrystalline  graphite  were  measured  in  the  dilatometer. 

Because  the  probes  are  also  made  of  poly crystalline  graphite,  the  net 
thermal  strain  measured  by  the  LVDT  was  predicted  by  Eq.  (1)  to  be  *  zero 
except  for  such  fine-scale  nonlinearities.  Uncorrected  [i.e.,  Eq.  (1)) 
expansion  curves  that  illustrate  this  behavior  were  obtained  from  three 
ATJ-S  graphite  samples  of  different  lengths  from  the  same  block  of  graphite 
and  are  plotted  in  Fig.  7.  Because  the  thermal  expansion  of  the  gruphite 
in  the  samples  and  probes  is  assumed  to  be  very  smooth,  the  inflected  shape 
of  the  curve  is  attributed  to  interaction  of  the  thermomechanical  and 
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Ftg,  6.  (a)  Initial  dilatometer  measurement,  following  calibration,  on 

NIST  sapphire  standard  reference  material  SRM  738  (solid  line). 
The  published  data  for  this  material  are  displayed  pointwise  as 
open  circles.  The  deviation  between  measured  and  published  data 
is  also  shown  (dashed  line).  The  indicated  error  is  in  the 
optical  calibration  data  (>  700°C). 

(b)  Two  successive  runs  on  SRM  738  after  subtracting  deviation 
data  (extrapolated  to  3000°C)  from  graphite-standard  calibration 
curve.  (See  also  Fig.  5.)  Published  data  for  SRM  738  are  also 
shown,  as  open  circles. 
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Fig.  7.  Thermal  expansion  data  for  «TJ-S  poly crystalline  graphite  samples 
of  different  lengths.  Since  the  ATJ-S  graphite  has  almost  the 
same  CTE  as  the  sample  tube/pushrod  assembly,  the  uncorrected  data 
have  a  small  net  slope.  The  inflections  present  in  the  curves  are 
attributed  to  the  instrumental  response  of  the  dilatometer,  since 
the  samples  and  the  sample  probes  undergo  no  phase  changes  during 
heating.  The  dilatometer 's  response  to  changes  in  sample  length 
is  nonlinear  and  not  monotonically  varying  with  respect  to  sample 
length.  The  length  of  sample  for  each  curve  is  indicated. 
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electronic  systems  in  the  dilatometer.  Although  it  varies  with  sample 
length  in  a  complex  manner,  the  fine-scale  instrumental  response  can  be 
considered  constant  if  the  sample  length  varies  only  by  the  amount  typi¬ 
cally  associated  with  the  thermal  strain  of  high-temperature  materials 
(<20  *  10~^  °C“1).  [For  very  large  strains,  such  as  those  associated  with 
py roly zing  carbon  precursor  materials  (see  Fig.  12),  fine-scale  nonlinear¬ 
ities  can  be  neglected.]  Such  nonlinearities  are  included  in  the  instru¬ 
mental  response  term  [Eq.  (3)]  by  using  a  higher-order  polynomial  function 
to  fit  the  data  prior  to  subtraction,  and  a  seventh-order  polynomial  fit 
was  found  to  be  adequate.  The  room- temperature  sample  length  also  must  be 
held  constant;  we  therefore  standardized  the  room- temperature  sample  length 
at  1  cm.  Based  on  the  data  shown  in  Figs.  5  and  6,  the  resolution  of  the 
system  is  better  than  ±5  urn,  giving  a  typical  sensitivity  of  about 
0.5  *  10*6  “C*1. 

The  above  rationale  and  data-correction  procedures  are  also  applied  to 
analyze  data  from  dilatometer  runs  with  bar  samples  in  the  flexure  mode, 
using  the  arrangement  shown  in  Fig.  4b. 
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IV.  RESULTS 


In  this  section  we  briefly  describe  the  results  from  early  runs  of  the 
dilatometer  to  illustrate  the  utility  and  flexibility  of  the  system.  The 
discussion  is  also  intended  to  inform  about  the  general  thermal  strain 
behavior  of  carbonaceous  materials. 

Uncorrected  data,  representing  the  net  thermal  strain  of  the  dilatom- 
eter-sample  system  [Eq.  (1)3,  are  presented  in  Fig.  8  for  non-heat-treated 
PG  and  the  PG  standard  sample.  The  slope  is  negative  because  the  sample 
data  are  uncorrected.  The  progressive  stabilization  of  the  PG  (which  is 
manufactured  at  about  2400°C)  and  the  expansion  of  the  PG  standard  (which 
had  been  pre-annealed  at  2750°C)  are  indicated.  The  magnitude  of  the 
sharp,  irreversible  increase  in  the  basal-plane  dimension  of  the  non-heat- 
treated  sample  diminishes  progressively  as  crystalline  defects  are  annealed 
out  and  the  degree  of  graphitic  order  of  the  sample  approaches  a  limiting 
value10-12  with  increasing  time  at  high  temperature.  The  CTE  of  the 
pyrolytic  graphite  in  the  basal  plane  changes  little  as  a  result  of  the 
high-temperature  heat  treatment,  indicating  good  stability  of  the  sample’s 
CTE  with  respect  to  use  of  the  standard  at  high  temperature. 

Typical  expansion  data  taken  from  various  carbon  materials  are  3hown 
in  Figs.  9  and  10.  The  true  (corrected)  data  taken  both  perpendicular  and 
parallel  to  the  basal  plane  of  a  sample  of  pyrolytic  graphite.  Fig.  9, 
exhibit  the  large  thermal  strain  anisotropy  characteristic  of  graphite,  as 
well  as  the  irreversible  strains  associated  with  further-graphitizing  heat 
treatment,  past  the  original  temperature  of  manufacture. 10  The  data  in 
Fig.  10  are  for  representative  samples  of  the  three  principal  weave  geom¬ 
etries  of  C-C  composites:  one-,  two-,  and  three-dimensional  weaves.  They 
are  very  similar  in  CTE  because  their  thermal  strains  are  dominated  by 
their  constituent  carbon  filaments,  which  vary  in  CTE  by  only  about  a 
factor  of  2  from  that  of  basal -plane  graphite.11 


23 


AL/Lq  (jim/cm) 


TEMPERATURE  (°C) 


Fig.  8.  Uncorrected  thermal  strain  data  from  pyrolytic  graphite.  Curves 
(a)  and  (b)  are  from  successive  runs  on  as-received  samples; 
curves  (c)  and  (d)  are  from  the  in-house-calibrated  standard.  The 
CTE  of  the  pyrolytic  graphite  (excluding  its  initial  high- tempera¬ 
ture  pyrolysis  expansion)  is  relatively  constant  even  after  the 
large  (irreversible)  strains  that  accompany  the  initial  heatings 
above  the  temperature  of  manufacture  (about  2400°C). 


AL/L0  (ftmtcm ) 


Fig.  9.  True  thermal  strain  response  of  pyrolytic  graphite  in  (a)  basal - 
plane-normal  direction  (c-axis)  for  the  second  run,  and  (b)  and 
(c),  parallel  to  the  basal  plane  (a-axis)  for  the  first  and  second 
runs,  respectively.  The  irreversible  shrinkages  and  expansions 
occurring  in  the  c-  and  a -axes  above  2500°C  are  associated  with 
further  graphiti2ation.y  Curve  (c)  indicates  that  the  sample 
becomes  stable  with  respect  to  thermal  strain  after  sufficient 
heat  treatment. 
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Fig.  10.  Thermal  strain  of  several  carbon-carbon  composite  configurations: 
(a)  a  Cartesian  223  weave  three-dimensional  composite  manufac¬ 
tured  by  General  Electric;  (b)  a  unidirectional  C-C  composite 
made  with  T300  fibers  and  pitch-derived  graphite  matrix;  and 
(c)  a  two-dimensional  composite  made  from  satin  weave  P100  fibers 
in  a  pitch-derived  graphitic  matrix. 
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As  previously  mentioned,  the  dilatooeter  can  also  be  run  in  the  three- 
point-bend  mode.  In  this  arrangement,  the  sample  is  in  the  form  of  a  bar 
0.8  x  2.4  x  5  mm  (length).  The  pushrod  is  placed  on  the  center  of  the  span 
to  record  the  inelastic  deflection  as  a  function  of  temperature.  Those 
elastic  deflections  that  occur  are  of  the  order  of  a  few  microns  and  are 
neglected.  This  type  of  test  is  somewhat  uncommon;  therefore,  it  is 
appropriate  to  indicate  the  general  type  of  data  yield. 

The  inelastic  strain  of  samples  of  pyrolyzing  PAA,  phenolic,  and 
polyfurfural  alcohol  during  heat  treatment  is  shown  in  Fig.  1 1 .  The  strain 
behavior  reveals  that  there  are  two  temperature  ranges  in  which  significant 
amounts  of  strain  energy  are  absorbed  by  the  samples.  The  first  occurs  at 
low  temperatures,  essentially  from  the  cure  temperature  to  about  1000°C, 
and  relates  to  the  formation  of  the  polymeric  carbon  submicrostructure . 1 ^ 
The  second  range  occurs  above  about  2000°C,  where  graphi tization  is  ini¬ 
tially  observed  in  the  (formerly  glassy)  matrices  of  C-C  composites,  up  to 
the  temperature  limit  of  the  dilatometer.  Macroscopic  work  absorption  in 
this  range  is  probably  associated  with  creep  and  grain-boundary  sliding 
processes. ^  Such  observations  are  important  in  research  on  carboniza¬ 
tion  and  graphitization  processes.  For  further  information,  the  reader  is 
referred  to  the  literature. 
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Fig.  11.  Data  from  three- point- bend  samples  of  (a)  P Aft,,  (b)  phenolic,  and 
(c)  furfural  resin  samples  that  were  precured  to  350°C.  The  max¬ 
imum  shear  and  outer-fiber  tensile  stresses  were  3.5  KPa  and  0.25 
MPa,  respectively.  The  data  indicate  that  the  three  carbons  dif¬ 
fer  in  the  amount  and  temperature  range  of  strain-energy  absorp¬ 
tion.  Phenolic  resin  absorbs  more  strain  energy  during  the  early 
stages  of  pyrolysis;  PkA  absorbs  more  at  graph it izat ion  tempera¬ 
tures  ( >2000*0);  and  furfural  resin  carbon  absorbs  less  at  all 
temperatures. 
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V.  CONCLUSIONS 


The  dilatometer  system  described  here  can  measure  the  thermally 
induced  strain  (deflection)  of  samples  to  a  precision  of  better  than  ±5  um 
from  room  temperature  to  3000°C.  This  precision  is  achieved  by  character¬ 
izing  the  instrumental  response  of  the  system  for  a  known  standard  that  can 
be  referenced  to  NBS  (NIST)  data,  then  algebraically  subtracting  that 
functionality  from  the  strain  of  the  unknown  sample.  The  dilatometer  thus 
emulates  a  virtual  two-pushrod  dilatometer  system.  Using  computer  data 
acquisition  and  analytical  methodology  makes  quantifying  a  sample's 
thermodynamic  responses  a  very  rapid  process. 

Non thermodynamic  properties,  such  as  the  strains  occurring  in  carbon¬ 
aceous  materials  during  their  pyrolysis,  can  also  be  measured.  Differences 
in  the  pyrolysis  behavior  of  carbon  precursors  have  been  observed.  It  is 
proposed  that  characterizing  and  analyzing  such  differences  may  shed  light 
on  the  behavior  of  precursors  in  the  matrix  of  C-C  composites  and  the 
effects  of  the  precursors  on  the  bulk  thermal  strains  of  such  composites. 
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APPENDIX  A. 


DILATOMETER  DATA  ACQUISITION  PROGRAM 


The  following  BASIC  language  code  runs  on  a  multitasking  Hewlett' 
Packard  personal  computer  under  a  BASIC  Shell  in  the  UNIX  environment.  It 
is  used  for  laboratory  data  acquisition. 
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10  . . mm... . M...MM 

20  !  •• 

30  !  ASTRO  Furnace  data  acquisition  program  *• 

40  !  •• 

50  !  Jun  16  1987  (V.rs  2.1)  ** 

80  !  •• 

70  !  •• 

80  !  Makes  a  daughter  file  in  /usr/ASTRO  •• 

90  !  which  can  be  read  by  Unix  for  VAX  transfer  •• 

100  !  •• 

110  !mm...mm.**..mmm..m.m..mmm.mmm*m..mm».m 
120  SHORT  ainT, maxT, inter*, tia«e,hrt,thick,sthick, aero  8  INTEGER  nopoints 
9  point=l 

130  CLEAR  8  ALPHA  0,0  9  DIM  saaplt[80] 

140  DISP  9  DISP  9  DISP  9  DISP 

150  DISP  •  Astro  Furnace  Data  program*  9  MASS  STORAGE  IS  */u»r/ASTR0* 
160  DISP  9  ASSIGN  16  TO  *astro! .dat.tap* 

170  DISP  •  calibrate  LVDT  [y/n] ■ 

180  INPUT  yesnot  9  PRINTER  IS  16 

190  IF  yesnof=»n*  OR  yesno*='N*  THEN  GOTO  220 

200  GOSUB  990 

210  GOTO  350 

240  GOTO  330 

330  DISP  9  DISP  9  DISP  •  enter  cal  factor  (Y/ua)* 

340  INPUT  calfae  0  sero=0 

350  DISP  9  DISP  8  DISP 

360  DISP  *  Enter  heating  rate  (C/ain.)1 

370  INPUT  hrt 

380  DISP  9  DISP  *  Enter  File  name* 

390  INPUT  naaeS 
400  sl=LEN (naaeS) 

410  dateeS=DATBS 

420  DISP  9  DISP  *  Saaple  I.D.  1 
430  INPUT  saaplt 

440  DISP  9  DISP  *  Enter  saaple  thickness  (aa)' 

450  INPUT  sthick 

460  DISP  9  DISP  1  Enter  ain.aax  temperatures  (C)* 

470  INPUT  ainT.aaxT 

480  DISP  9  DISP  *  Enter  saaple  interval  (seconds)* 

490  INPUT  interr 

500  DISP  9  DISP  9  DISP  *  Length  of  run  (minutes)1 
510  INPUT  tiaee 

520  tiaeee=(aaxT-aisT)/hrt  !  time  length  of  run  !!!! 

530  DISP  9  DISP  *  Time  to  aaxiaua  temp:  *;tiaeee  9  DSP 

540  nopnts=INT (tiaee/ (interr/60))  !  total  #  of  data  points  in  run  ! 

560  GOSUB  890  !***********************************8tart  file************ 
570  DISP  9  DISP  0  DISP  8  DISP  *  to  start  run,  Press  fl* 

580  ON  KEY#  1,*  start*  GOTO  600 
590  GOTO  590 

600  OFF  KEY#  1  9  DISP  9  DISP  'starting*  9  DISP  9  DISP 
9  DISP  *  Point  tiae  teap  LVDT* 

610  ON  KEY#  2,*  HOLD*  GOTO  1380 

620  tima=TIME  !  starting  seconds  since  midnight  ! 

630  datt=DATB- 87 000  !  starting  #  of  days  since  start  of  year 
640  !!!!!!start  of  data  loop! I 
650  ! 

660  ON  ERROR  GOTO  690 
670  ASSIGN  7  TO  *hpib.b2* 

680  OFF  ERROR  9  GOTO  695 
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690  WAIT  900  •  WTO  670 
695  OK  BU01  GOTO  recorer 

700  OUTPUT  709  ;  aAC5a  9  KKTSK  709  ;  L1,L2,L3,U,L5,L6 
710  OUTPUT  709  ;  aAC8a  9  EWER  709  ;  T1,T2,T3,T4,TS,T6 
720  OFF  ERROR  •  ASSIGN  7  TO  ••• 

730  ! 

740  lTdt*(Ll*L2*L3*L4-*-L5-*-L6) /8/c*lf *c  I  t**p-(Tl*T2*T3*T4*T&*T8)/6*1000 

760  IF  DATt-S7000xUtt  TIEN  elti*s-(TDIE-ti**)/60  «  GOTO  780 

780  IF  DATE- 87000- datt-1  TIBI  slti**«(TIMR*ti**)/eo  •  GOTO  780 

770  IF  DAIE-87000-d*tt>l  TIBI  eltiM-(TIMR*tia*)  /60*  (DATE-87000-datt)  .1440 

780  C-CBIL (point* 15) 

790  i 
800  1 

810  DISP  point;*  a;INT(*lti***1000)/1000;a  MNT(t«*f>*lOO)/100; •  •; 

INT  (lvctt*  100)  /100 

820  point*point*l  8  IF  point>nopata  TIBI  WTO  860 

830  PUNT  point ;  BIT (Irdt* 1000) /1000 ,INT(t*ap* 1000) /10OO,INT(#lti»#* 1000) /1QOO 
836  ! 

840  HIT  intorr *1000-5000 
850  GOTO  650  !  n«xt  point  t 

860  nsa*8[sl+l]«a.dsta  9  RENAME  a**tro! .dnt.tap*  TO  naael 
870  DISP  0  DISP  0  DISP  *  BID...  (bye)* 

880  END 

890  !  SIAM'S  FILE  . . . . . 

900  !  •• 

960  PUNT  naaeS,d»t**S,*inT,*»xT  9  PUNT  ssaplt  I  •• 

970  PUNT  hrt,sthick,sero,TIMB  9  PUNT  'point#  *;*  Micron*  a; 

*  Teaperature  elapsed  tiae* 

980  KETUKN 
990  !  calibrates  LVDT 
1000  ! 

1010  ! 

1020  ! 

1030  DISP  9  DISP  9  DISP  *  Z«ro  LVDT,  then  push  rot  urn. ! 

1040  INPUT  duaay*  9  DISP  9  DISP  1 
1050  ON  KUOK  GOTO  1080  ! 

1060  ASSIGN  7  TO  ahpib.b2a  t 
1070  GOTO  1090  ! 

1080  WAIT  200  9  GOTO  1060  ! 

1090  OFF  EE&OR  ! 

1100  FOR  LL=1  TO  3  ! 

1110  OUTPUT  709  j  aAC5a \ 

1120  ENTER  709  ;  *1, *2, *3, *4, *5, *6, *7, *8, *9  ! 

1130  VAR  1000  ! 

1140  eeroo (LL)=(sl*s2**3*s4+*5*»6*e7fs8*i9) /W  ! 

1150  NER  LL  ! 

1160  s*ro»(Mro(l)+s«ro(2)+**ro(3))/3  ! 

1170  DISP  9  DISP  *  pine*  shia  in  fixture,  *nt*r  it*  tbickn****  ! 

1180  DISP  *  and  pres*  return1  ! 

1190  DISP  9  DISP  a  Shi*  tbickne**  (*icron*)a  I 
1200  INPUT  thick  I 

1210  DISP  9  DISP  9  DISP  a  (...press  RETURN  when  tbe  shin  is  in)1  ! 

1220  INPUT  doasyf  ! 

1230  ON  ERROR  GOTO  1260  ! 

1240  ASSIGN  7  TO  ahpib.b2a  ! 

1250  GOTO  1270  I 
1280  VAR  200  9  GOTO  1240  ! 

1270  OFF  ERROR  ! 

1280  FOR  LUl  TO  3  I 
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1200  OUTPUT  700  ;  »AC5*  ! 

1300  ENTS*  700  ;  cl,c2,c3)c4,c5,c6,c7Jc«lc0  ! 

1310  e*lis(LL)*((cl+c2*c3+c4+c£'»'c6'»e7*c8+c0)/0-s«ro)/tliidi[  !  Y/ua 
1315  »AIT  1000  ! 

1320  NEXT  LL  ! 

1330  calfac*(ealfa(l)+calfa(2)+calfa(3))/3  ! 

1340  ASSIGN  7  TO  •••  ! 

1350  DISP  0  DISF  *  Cal.  factor** ; calf ac  0  RETURN  ! 

1360  ! - 

1370  ! - 

1380  DISP  8  DISP  8  DISP  *  ...holding,  #nd*i3;  coatinue*f4* 
1300  ON  KEY#  3,»  end*  GOTO  860  8  ON  KEY#  4,*  coat.*  GOTO  1410 
1400  GOTO  1400 

1410  OFF  KEY#  3  8  OFF  KEY#  4  8  GOTO  840 
1420  recover: 

1430  DISP  *  A  Tineout  has  Occurred  on  701 . . .  Recovering* 

1440  ASSIGN  7  TO 
1450  GOTO  650 


uu 

WWFW 

uu 

Iff  I 
#### 
#### 
#### 

nil 

wTVW 

Iff# 

ffff 

Hff 

ffff 
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APPENDIX  B. 


DILATOMETER  DATA  ANALYSIS  PROGRAM 


The  following  IDL  language  code  i:  -un  on  a  VAX  11/785  computer  to 
analyze  the  data  collected  in  the  laboratory.  Data  are  transferred  to  the 
VAX  via  Kermit  file-transfer  software  (public  domain). 


^Interactive  Data  Language,  RSI,  Inc.,  2021  Albion  St.,  Denver,  CO. 
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SET  PLOT,0 

VTIOO 

!FANCY=0 


Version  S.OOpyrostd  baseline  7th  polynomial 
JAN  10,  1989  *ith  SMOOTHING 
and  SCATTER  PLOT  and  AOTOREAD 
and  DERIVATIVE  and  LBGANDS 


PVAR=” 

print, ’  ASTRO  PILE  DATA  PLOTTER;  OVERPLOTTED  PILES  WILL  BE  OYERLAYED  IN  THE1 
PRINT, ’  ORDER  GIVEN.  SCALE  WILL  BE  THAT  OF  THE  P3SST  PILE  GIVEN. ’ 

PRINT,’  Be  sure  that  the  terminal  is  in  CAPITAL  letter  mode! 

PRINT,’  ’ 

PRINT,’  with  SMOOTHING  and  SCATTER  PLOT’ 

PRINT,’  ’ 

NUMP=10 

PRINT,’  Enter  #  o f  points  to  be  averaged 
read.NUMP 
PRINT,’  ’ 

PRINT, ’  NUMBER  OP  FILES  TO  BE  PLOTTED  ? 

AB=2 
Z=1 

RBAD.AB 
C0DE=1 


for  smoothing’ 


(<9)  ’ 


;  NUMBER  OF  PILES 
;  A  FLAG 


0RD=2 
PRINT,’  ’ 

PRINT, ’  Plot  a  polynomial  fit  to  the  data?  [Y/N] ’ 
read , PVAR 
pjjjjT  y  * 

IF  PVAR  EQ  ’Y’  THEN  BEGIN 
PRINT,’  WHAT  POLYNOMIAL  ORDER  ?’ 

READ,  (3RD 
ENDIF 


PRINT,’ 

PRINT,’ 

PRINT,’ 

> 

X-AXIS 

Y-AXIS 

CODB  ’ 

1 

PRINT, ’“ 
print, ’ 
PRINT,’ 

Temperature 

Time 

Time 

Expansion 

Expansion 

Temperature 

1’ 

2’ 

3’ 

PRINT,’ 

Temperature 

Derivative  (Temp.) 

4 

(for 

dilation) ’ 

PRINT,’ 

Time 

Derivative  (Time) 

5 

(for 

dilation) ’ 

PRINT,’ 

Temperature 

Deflection 

6 

(for 

flexure) ’ 

PRINT,’ 

Temperature 

Deflection 

7 

(for 

compression) ’ 

PRINT,’ 

Temperature 

Derivative 

8 

(for 

compression) ’ 

PRINT,’ 

C0DR=2 

Temperature 

Derivative 

8 

(for 

flexure) ’ 

PRINT,’  ’ 

PRINT, ’  Enter  CODE’ 
print, ’  ’ 

READ, CODE 


NUM=INDGEN(AB) 

NAME*STRARR(20,AB) 

NAMEE=” 


;  VECTOR  OF  POINT  NUMBERS 
CORRESPONDING  VECTOR  OF  NAMES 


YSNO«” 

READ, ’Display  the  Directory  ' ,YSNQ 
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IP  YSNO  EQ  ’Y’  TEEN  SPAIN,  ’DU’ 

NUMBERS 

FOB  1*1, AB  DO  BEGIN  ;  get  oases  of  files  to  be  plotted 

PRINT,  ’ENTER  NAME  OP  PUB  ’ 

READ, NAMES 
NAME (1-1) *NAMEE 
endfor 

STD1*  [-4 . 35479E-15, 4 . 77169E-12 ,  -1 . 9023SB-9 , 3 . 18491E-7 ,  -1 . 44959E-5 ,  -1 . 18590B-3  ,-.44418,1.34 

;  BASELINE  PIT  DATA  JUN.  25,  1988 
STD2* [-3 . 95873E-16, 3 . 98258B-12 , -1 . 56485B-8 , 3 . 14101B-5 , -1 . 35172B-2 , . 31399] 

;  COEPF’S  FOR  GRAPHITE  IN  BASAL  PLANE 
STD3* [6 . 66491E-15, -5 .51405B-11 , 1 . 69794B-7 , -2 . 40729E-4.0. 13018] 

;  BASELINE  FOR  COMPRESSION  SAMPLE 

for  i*l,ab  do  begin  ;  do  plotting... 

0PENR,1,NAME(I-1) 

A*STRARR(30,1) 

B=STRARR(60,1) 

C*flt»rr(4) 

D=” 

E=DOUBLE (FLTARR (4 , 2000) )  ;  MAKE  DATA  VECTOR  DOUBLE  PRECISION 

ON  IOBRROR, POINTS 

READF  l  A  B  C  D  E 

POINTS: 'PRINT) (MIN(WHERE(E  EQ  0)))/4,’  POINTS  BEAD’ 

ON  IOERROR,NULL 

PRINT, ’  Nuaber  of  Data  Points  to  be  Plotted  ? ’ 

READ,NUMEE 

NUM(I-1)=NUMEE 

close, 1 


X=E (2,0: NUM (1-1) ) 
Y=E(1,0:NUM(I-1)) 
XT=E(3,0:NUM(I-1)) 


DEFINE  NEW  VARIABLES  AND  SET  UP  PLOT 

MICRONS 

TIME 


PYROFLX* (-6 . 6E-3) *X 


FOR  FLEX  SAMPLES 


Y*Y-  (Y  (0)  +Y  (1 )  +Y  (2)  )  /3 


SET  Y(B.T.)=0 


X*X/8 

Y1=STD1 (0) *X*7+STD1 (1) *X*8+STD1 (2) *I*5+STD1 (3) *I*4+STD1 (4) *X‘3>STDl (5) *X‘2*STD1 (8) *X+STD 

;  MAKE  BASELINE  DATA 

X*X*8 

Y2=STD2  (0)  *X“5+STD2(1)  *X‘4+STD2  (2)  *X“3+STD2  (3)  *X*2+STD2  (4)  *I*STD2(5) 

•  MAKE  BASAL  PLANE  DATA 
;  FOR  GB  STANDARD 


IF  (CODE  SQ  7)  OR  (CODE  EQ  8)  THEN  BEGIN 

Y4-STD3 (0) *X“5*STD3 (1) *X‘4+STD3 (2) .X‘3*STD3 (3) *X*2*STD3 (4) *X 

;  MAKE  BASELINE  FOR  COMPRESSION 

END  IF 


H2 


!GRXD*2 

!XTITLE* ’Teaperature  (C)  ’ 
lYTITLE* ’Delta  L/Lo  (aicrons/ca) ’ 


!MTITLE=’ THERMAL  EXPANSION  Be  atmosphere’ 
IF  CODE  EQ  2  THEN  !XTITLE=’TIME  (minutes) * 
IF  CODE  EQ  3  THEN  BEGIN 
!XTITLE=’TIME  (minutes) * 

! YTITLE* * TEMPERATURE  (C) ’ 

ENDIF 

IF  CODE  EQ  4  THEN  BEGIN 

!  YTITLB=  ’dS/dT  (ua/mn/deg .  C)  ’ 

!MTITLE= ’HIGH  TEMPERATURE  FLEXURE’ 

ENDIF 

IF  CODE  EQ  9  THEN  BEGIN 

!MTITLE=’HIGH  TEMPERATURE  FLEXURE’ 
!YTITLE=’ dS/dT  (um/mm/C) ’ 

ENDIF 

IF  CODE  EQ  8  THEN  BEGIN 

!MTITLE=’HIGH  TEMPERATURE  COMPRESSION’ 
!YTITLE=*dL/dT  (u*/aa/C)  ’ 

ENDIF 

IF  CODE  EQ  5  THEN  BEGIN 

!YTITLE=’dS/dt  (ua/ca/min.) * 

! XTITLE= ' Time  (minutes) ’ 

!MTITLE=’HIGH  TEMPERATURE  FLEXURE’ 

ENDIF 

IF  CODE  EQ  6  THEN  BEGIN 
!YTITLE=’ Deflection  (urn)’ 

! XTITLE= ’ Temperature  (C) ’ 

!MTITL£=’HIGH  TEMPERATURE  FLEXURE’ 

ENDIF 

IF  CODE  EQ  7  THEN  BEGIN 

!MTITLE=’HIGH  TEMPERATURE  COMPRESSION’ 

! YTITLE= ’ Compressive  Strain  (10*-3) ’ 
ENDIF 


IF  (CODE  EQ  1)  or  (code  eq  4)  THEN  BEGIN 
Y=(Y-Y1+Y2)/(C(1)*10)  ;  CORRECT  DATA  FOR  DILATION  SAMPLES 

ENDIF 


IF  (CODE  EQ  7)  OR 
Y=(Y-Y4)/C(1) 

ENDIF 


(CODE  BQ  8)  THEN  BEGIN 

j  SUBTRACT  COMPRESSION  BASELINE 
;  AND  MAKE  INTO  um/mm 


IF  (CODE  BQ  6)  OR  (CODE  EQ  9)  THEN  BEGIN 
Y=Y-PYROFLX  ;  SUBTRACT  FLEXURE  BASELINE 

ENDIF 


SBTVIEHPORT,  .11,  .80,  .22,  .9 


IF  Z  GT  1  THEN  BEGIN  ;  USE  ’PLOT’ 

!LINETYPE=I-1 

IF  (CODE  EQ  1)  or  (CODE  EQ  8)  OR  (CODE  EQ  7) 

IF  CODE  EQ  2  THEN  OPLOT, XT, SMOOTH (Y.NUMP) 

IF  CODE  BQ  3  THEN  OFLQT,XT,X 
IF  (CODE  EQ  4)  OR  (CODE  EQ  8)  OR  (CODE  EQ  9) 
Dy=DERIV (X, (SMOOTH (Y,NUMP) ) ) 

Dy(0)*Dy(3) 

Dy(l)=D7(0) 

OPLQT,X,Dy 


ONLY  IF  1ST  TIME  AROUND 

THEN  QFLOT,X, SMOOTH (Y,NUMP) 
;  ELSE  USE  ’OPLOT’ 

;  TIME: EXPANSION 

THEN  BEGIN 


END  IF 

IF  CODS  BQ  5  THEN  BEGIN 
Dy*DEiIY (IT ,  (SMOOTH(Y,NUMP))) 

Dt(0)«Dt(3) 

Dr(l)-l>y(0) 

OPLOT, XT, Dr 
END  IF 

ENDIF  ELSE  BEGIN 

!PSYM=3  ;  WJ  DOTS 

VTIOO 

1LINETYFB-0 

PL0T,X,Y 

IF  (CODE  EQ  1)  OR  (CODE  EQ  6)  OK  (CODE  BQ  7)  TEEN  PLOT,X(Y 

IF  CODE  BQ  2  THEN  PLOT, XT, Y 

IF  CODE  BQ  3  THEN  PLOT, XT, X 

IF  (CODE  EQ  4)  OR  (CODE  BQ  8)  OR  (CODE  BQ  6)  THEN  BEGIN 
Dp4)ERIY (X,  (SMOOTH (Y , NUMP) ) ) 

Dy(0)*Dy(3)  ;  GENERAL  DERTVATTYB  CURVES 

Dr(l)*Dy(0) 

PLOT, X, Dr 

ENDIF 

IF  CODE  EQ  5  THEN  BEGIN 

Dr=DERIV (IT , (SMOOTH (Y.NUMP)}) 

Dr (0)*Dy (3) 

Dr(l)-Dy(0) 

PLOT, XT, Dr 

ENDIF 
!PSYM=0 

IF  PVAR  NE  »Y»  THEN  BEGIN 
SET  XY 
!PS?M=0 
!LINETYPE=0 

IF  (CODE  EQ  1)  OR  (CODE  BQ  6)  OR  (CODB  BQ  7) 

IF  CODE  EQ  2  THEN  OPLOT, XT, SMOOTH (Y,WJMP) 

IF  CODE  EQ  3  THEN  OPLOT, XT, I 
IF  (CODE  EQ  4)  OR  (CODE  EQ  8)  OR  (CODB  BQ  0) 

Dr*DERIV(X, (SMOOTH (Y.NUMP))) 

Dy(0)»Dr(3) 

Dr(i)*Dr(0) 

OPLOT, X, Dr 
ENDIF 

IF  CODE  BQ  5  THEN  BEGIN 

Dr*DERIV (XT, (SMOOTH (Y.NUMP))) 

Dr(0)-Dr(3) 

Dr(i)*Dr(0) 

OPLOT, XT, Dr 
ENDIF 

ENDIF 

IF  PVAR  BQ  *Y’  THEN  BEGIN 
COEFF*POLY  FIT(X,Y,ORD,GGG) 

IPSYM-0  “ 

OPLOT,!, GGG 
ENDIF 


; OPLOT  SMOOTHED  CURVE 

THEN  OPLOT,!, SMOOTH (Y.NUMP) 
THEN  BEGIN 


;  OPLOT  POLYNOMIAL  CURVE  IF 
;  PVAR-Y 


Z-2 

BNDELSE 


BNDFOR 

!Noms»i 

SET  VIBWPORT.O, 1,0,1 
SET“XY,0, 1,0,1 
ABC5=.82 

XYOUTS,  .85,  .86,  ’Hums:  ' 

FOE  LL=1,AB  DO  BEGIN 
NNN=FINDGEN(AB*1) 

!LINETYPE=NNN(LL-1) 

OFLOT, [.82, .88] , [ABCD-(LL-1)*.05,ABCD-(LL-1)*.05] 
XYOUTS,. 875, ABCD-(LL-1)*.05-. 01, NAME (LL-1) 

ENDFOB 
SET  XY 
SET" VIEWPORT 
!NOERAS=0 

SET  VIEWPORT, 0, 1,0,1 
SET  XY, 0,1, 0,1 

XYOtJTS,  .17,  .1,  ’Nue,  Run  Date  :  ’ 

XYOUTS, .30, .1, A 

XYOUTS, . 17, .05, ’Heating  Rate  (C/min.) :  ’ 

XYOUTS, .455, .05, C(0) 

XYOUTS, .7, .05, ’Lo  (ca) :  * 

XYOUTS, .79, .05,0(1) 

XYOUTS,. 17,. 075, ’Saaple  ID:  ’ 

XYOUTS, .37, .075, B 

SBT_XY 

SETVIEWPORT,  .11,  .80,  .22,  .9 


;  ****** ************************** a************** ************ »•*****• 

; PLOTTING  SECTION 

;  **********************************************  ***•***••••******••** 


PRINT,  ’  SAVE  THIS  ?  [Y/N] » 

YN=” 

BEADjYN 

IF  YN  EQ  ’Y’  THEN  BEGIN  ;  IF  HARD  COPY  IS  DESIRED 

Z=1 

for  isl,ab  do  begin  ;  do  plotting... 

OPENR, 1 ,NAMB(I-1) 

A=STRARR(30,1) 

B=STRARR(60,1) 

C=f ltarr (4) 

D=” 

E=D0UBLE(FLTABR(4,NUM(I-1)+1))  ;  MAKE  DATA  VECTOR  DOUBLE  PRECISION 

READF,1,A,B,C,D,E 

CLOSE, 1 


X=E(2,0:NUM(I-1)) 


;  DEFINE  NEW  VARIABLES  AND  SET  UP  PLOT 


Y«B(1,0:NUM(I-1))  ;  MICRONS 

XT*B(3,0:NUM(I-1))  ;  TIME 

PYROFLX=  (-6 . 8E-3)  *1 

;  FOR  FLEX  SAMPLES 

Y*Y- (Y(O)  +Y  (1) +Y (2) ) /3  ;  SET  ZERO  OF  CURVE 

X=X/8 

Y1=STD1 (0) *X*7+STD1 (1) *X‘6+STD1 (2) •X*5+STDl (3) *I*4+STD1 (4) *X‘3+STD1 (5) *X‘2-STT>1 (8) .X+STD1 ( 

;  MAKE  BASELINE  DATA 

X=X*8 

Y2=STD2 (0) *X*5+STD2 (1) *I*4+STD2 (2) *X‘3+STD2 (3) *X‘2+STD2 (4) *X*STD2 (5) 

;  MAKE  BASAL  PLANE  DATA 
;  FOR  GR  STANDARD 


IF  (CODE  BQ  1)  or  (code  eq  4)  THEN  BEGIN 
Y=(Y-Y1+Y2)/(C(1)*10)  ;  CORRECT  DATA  FOR  DILATION  SAMPLES 

ENDIF 


IF  (CODE  EQ  7)  OR  (CODE  EQ  8)  THEN  BEGIN 
Y=(Y-Y4)/C(1)  ;  SUBTRACT  COMPRESSION  BASELINE 

;  AND  MAKE  INTO  ua/aa 


ENDIF 


IF  (CODE  EQ  6)  OR  (CODE  EQ  0}  THEN  BEGIN 
Y=Y-PYROFLX  ;  SUBTRACT  FLEXURE  BASELINE 

ENDIF 


.'fancy  =  2 
set_j>lot,4 

SET_VIEWP0RT,.11,.80,.22,.9 

IF  Z  CT  1  THEN  BEGIN  ;  USE  ’PLOT’  ONLY  IF  1ST  TIME  AROUND 

j linetype=i~i 

ip  (CODE  EQ  1)  or  (CODE  EQ  6)  OR  (CODE  BQ  7)  THEN  OPLOT.X, SMOOTH (Y.NUMP) 

;  ELSE  USB  ’OPLOT’ 

IF  CODE  EQ  2  THEN  OPLOT, XT, SMOOTH (Y,NUMP)  ;  TIME: EXPANSION 

IF  CODE  EQ  3  THEN  OPLOT, XT, X 

IF  (CODE  EQ  4)  OR  (CODE  EQ  8)  OR  (CODE  EQ  9)  THEN  BEGIN 
Dy=DERIV (X, (SMOOTH (Y,NUMP))) 

Dy(0)=Djr(3) 

Dy(l)=Djr(0) 

OPLOT, X,Dy 
ENDIF 

IF  CODE  EQ  S  THEN  BEGIN 
Dy=DERIV(XT1 (SMOOTH (Y.NUMP))) 

Dy(0)=Dy(3) 

Dy(l)=Dy(0) 

OPLOT,XT,Dy 

ENDIF 

ENDIF  ELSE  BEGIN 

!PSYM*3  ;  DO  DOTS 

VT100 

!LINETYPE=0 

PLOT,X,Y 

IF  (CODE  EQ  1)  OR  (CODE  BQ  6)  OR  (CODE  BQ  7)  THEN  FLOT,X,Y 
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IF  CODE  EQ  2  THEN  PLOT, IT, Y 

IF  CODE  EQ  3  THEN  PLOT, XT, X 

IF  (CODE  EQ  4)  OK  (CODE  EQ  8)  OB  (CODE  EQ  9}  THEN  BEGIN 
Dy=DERIV(X, (SMOOTH(Y,NUMP))) 

Dy(0)=Dy(3)  ;  GENERAL  DERIVATIVE  CURVES 

Dy(l)-Dy  (0) 

PLOT,X,Dy 

ENDIF 

IF  CODE  EQ  5  THEN  BEGIN 

Dy=DERIV (XT , (SMOOTH (Y.NUMP))) 

Dy(0)=Dy(3) 

Dy(l)=Dy(0) 

PLOT, XT, Dy 

ENDIF 

!PSYM=0 

IF  PVAR  NE  *Y’  THEN  BEGIN 
SET  XY 
!PS?M=0 
!LINETYPE=0 

IF  (CODE  EQ  1)  OB  (CODE  EQ  8)  OB  (CODE  BQ  7) 

IF  CODE  EQ  2  THEN  OPLOT, XT, SMOOTH (Y,NUMP) 

IF  CODE  EQ  3  THEN  OPLOT, XT, X 
IF  (CODE  EQ  4)  OB  (CODE  EQ  8)  OB  (CODE  EQ  9) 

Dy=DERIV(X, (SMOOTH (Y.NUMP))) 

Dy(0)=Dy (3) 

Dy(l)=Dy(0) 

OPLOT, X,Dy 
ENDIF 

IF  CODE  EQ  5  THEN  BEGIN 

Dy=DERIY(XT, (SMOOTH (Y.NUMP))) 

Dy(0)=Dy(3) 

Dy (l)=Dy(0) 

OPLOT, XT, Dy 
ENDIF 

ENDIF 

nOUTS, !CXMAX*0.843, !CYMAX*1.02, ISTIMB 

IF  PVAR  EQ  ’Y’  THEN  BEGIN  ;  OPLOT  POLYNOMIAL  CURVE  IF 

COEFF=POLY_PIT (X , Y , OBD , GGG)  ;  PVAR=Y 

!PSYM=0 
OPLOT, X, GGG 

ENDIF 


; OPLOT  SMOOTHED  CURVE 

HEN  OPLOT, X, SMOOTH (Y.NUMP) 
THEN  BEGIN 


Z=2 

ENDELSE 

ENDFOB 

endif 

IF  CODE  EQ  1  THEN  BEGIN 
IF  Z  EQ  2  THEN  BEGIN 
IF  AB  EQ  1  THEN  BEGIN 

PRINT, *  TO  CALCULATE  ALPHA,  ENTER  NUMBER  OF  REGIONS  TO  BE  FITTED. ’ 
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PRINT,’  TO  END,  ENTER  ZERO’ 

ANSWE&=0 
READ, ANSWER 

if  ANSWER  GT  0  THEN  BEGIN 
FOR  1*1, ANSWER  DO  BEGIN 
PRINT,’  ’ 

PRINT, ’  ENTER  TEMP  RANGE  (LOWER, UPPER) ’ 

L0WER*5 

UPPER*5 

READ, LOWER, UPPER 

LOWER=MIN (WHERE (X  GT  LOWER))  ;  GET  SUBSCRIPTS  OF  TEMP  RANGE 

UPPER=MAX (WHERE (X  LT  UPPER)) 

XFIT=X (LOWER :UPPEfi)  ;  DEFINE  NEW  DATA  VECTORS  FOR 

YFIT*Y (LOWER : UPPER)  ;  ALPHA  FIT 

COEFF=POLY  FIT (XFIT,YFIT,1, LINE)  j  FIT  LINE  TO  DATA 

!LBQSTYPE=8 

OPLOT,XFIT,LINE 

ALPHA=COEFF ( 1 ) ♦ 100  ;  MAKE  ALPHA  INTO  MICRONS  PER  METER 

XYOUTS,X (UPPER) ,Y (UPPER) , ALPHA 

SETPLOT.O 

SET_VTBWPORT,  .11,  .80,  .22,  .9 

!FANCY=0 

OPLOT , XFIT , LINE 

XYOUTS.X (UPPER) ,Y (UPPER) .ALPHA 

SET_PL0T,4 

SET_V3EWP0RT, .11, .80, .22, .9 

!FANCY=3 

BNDFOR 
ENDIF 
END  IF 

ENDIF 

ENDIF 

IF  YN  EQ  *Y*  THEN  BEGIN 


!N0BRAS=1 

SET  VIEWPORT, 0,1, 0,1 

SET"*XY,0,1,0,1 

ABCfc=.82 


XYOUTS, .85, .86, ’Naaes: ’ 

FOR  LL*1,AB  DO  BEGIN 
NNN=FINDGBN (AB+1) 

!LINETYPB=NNN (LL-1) 

OPLOT, [.80, .87] ,  [ABCD- (LL-1) * . 05 , ABCD- (LL-1) * . 05] 
XYOUTS, .875, ABCD- (LL-1) *.05-. 01, NAME (LL-1) 

BNDFOR 


XYOUTS,. 17,. 1,’Naae,  Run  Date:  ’ 

XYOUTS,. 31,. 1, A 

XYOUTS , . 17 , . 05 , ’Heating  Rate  (C/ain . ) :  ’ 
XYOUTS, .33, .05,0(0) 

XYOUTS,. 6,. 05, ’Lo  (ca):  ’ 


XY0UTS, .87, .05,0(1} 

XYOUTS, .17, .075, ’S*apl«  ID:  * 

IYOUTS, .28, .075,8 

SET  XY 

SBOlBWPOET,  .11,  .80,  .22,  .8 
!NOSEaS-0 

YTlOO 

FEINT,  *  LASKB  PUNTING  THE  FILE  NOW. . ,  ’ 
WAIT, 2 
SET  PLOT,0 
VTl50 

SPAWN, ’LASE  DOWN  QMSPLOT . LIS ’ 

BNDIP 

SETPLOT.O 
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TECHNOLOGY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  "architect-engineer"  for  national  security 
programs,  specializing  in  advanced  military  space  systems.  The  Corporation’s  Technology 
Operations  supports  the  effective  and  timely  development  and  operation  of  national  security 
systems  through  scientific  research  and  the  application  of  advanced  technology.  Vital  to  the 
success  of  the  Corporation  is  the  technical  staffs  wide-ranging  expertise  and  its  ability  to  stay 
abreast  of  new  technological  developments  and  program  support  issues  associated  with  rapidly 
evolving  space  systems.  Contributing  capabilities  are  provided  by  these  individual  Technology 
Centers: 

Electronics  Technology  Center:  Microelectronics,  solid-state  device  physics, 

VLSI  reliability,  compound  semiconductors,  radiation  hardening,  data  storage 
technologies,  infrared  detector  devices  and  testing;  electro-optics,  quantum  electronics, 
solid-state  lasers,  optical  propagation  and  communications;  cw  and  pulsed  chemical 
laser  development,  optical  resonators,  beam  control,  atmospheric  propagation,  and 
laser  effects  and  countermeasures;  atomic  frequency  standards,  applied  laser 
spectroscopy,  laser  chemistry,  laser  optoelectronics,  phase  conjugation  and  coherent 
imaging,  solar  cell  physics,  battery  electrochemistry,  battery  testing  and  evaluation. 

Mechanics  and  Materials  Technology  Center:  Evaluation  and 
characterization  of  new  maceriais:  metals,  alloys,  ceramics,  polymers  and  their 
composites,  and  new  forms  of  carbon;  development  and  analysis  of  thin  films  and 
deposition  techniques;  nondestructive  evaluation,  component  failure  analysis  and 
reliability;  fracture  mechanics  and  stress  corrosion;  development  and  evaluation  of 
hardened  components;  analysis  and  evaluation  of  materials  at  cryogenic  and  elevated 
temperatures;  launch  vehicle  and  reentry  fluid  mechanics,  heat  transfer  and  flight 
dynamics;  chemical  and  electric  propulsion;  spacecraft  structural  mechanics, 
spacecraft  survivability  and  vulnerability  assessment;  contamination,  thermal  and 
structural  control;  high  temperature  thermomecbanics,  gas  kinetics  and  radiation; 
lubrication  and  surface  phenomena. 

Space  and  Environment  Technology  Center:  Magnetospheric,  auroral  and 
cosmic  ray  physics,  wave-particle  interactions,  magnetospheric  plasma  waves; 
atmospheric  and  ionospheric  physics,  density  and  composition  of  the  upper 
atmosphere,  remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared 
astronomy,  infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere;  effects 
of  electromagnetic  and  particulate  radiations  on  space  systems;  space  instrumentation; 
propellant  chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
atmospheric  chemical  reactions,  atmospheric  optics,  light  scattering,  stale-specific 
chemical  reactions  and  radiative  signatures  of  missile  plumes,  and  sensor  out-of-field- 
of-view  rejection. 


